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ABSTRACT 

Microquasars have been expected to emit high energy 7-rays due to their general sim- 
ilarities to the 7-ray emitting blazars (evidences of relativistic jets, non-thermal radio 
to X-ray emission). In fact, the first source of this type, Cyg X-3, has been recently 
unambiguously discovered by the satellite telescopes. We study the features of the 
7-ray radiation produced in these sources by relativistic electrons, accelerated in the 
inner part of the jet. The electrons initiate an Inverse Compton pair cascade in 
the radiation field of the accretion disk. Due to the anisotropy of the accretion disk 
radiation field, the spectra of 7-rays show strong dependence on the observation angle, 
the location of the emission region within the jet and the details of the acceleration 
process. As an example, we confront our model with the observations of the micro- 
quasar Cyg X-3, which has been recently reported as a transient GeV 7-ray source by 
the Agile and the Fermi Observatories. Satisfactory description of the 7-ray spectra 
observed from Cyg X-3 are obtained in the case of the injection of electrons in the 
inner part of the jet (located within 300 inner disk radius from the jet base) provided 
that the observer is located at relatively small angle to the jet axis. 

Key words: binaries: close: individual (Cyg X-3) — radiation mechanisms: non- 
thermal — gamma-rays: theory 



1 INTRODUCTION 

Microquasars are binary systems which show evidences of 
relativistic collimated outflows resembling those observed 



in active galactic nuclei (e.g. Mirabel fc R odriguez 1 1994 



iTingav et all Il995l : iHiellming fc Rupenl 1 19951 ). They are 



interpreted in the model of a jet produced in the in- 
ner part of the accretion disk around solar mass black 
hole whi ch accretes the matter from a close companion 
star ( e. g. iLevinson fc Blandfo rd 1996; Atovan fc Aharonianl 
19991 : iGeorganopoulos et al.l "|2002|: [Romero et all |2002|; 



Bosch-Ramon fc Par cdcs 20(3). Up to now, only one source 



of this type, i.e. Cyg X-3, has been clearly detected in 
the GeV 7-rays by the Agile and th e Fermi observatories 
|Tavani et aljfeooalAbdo et al.ll2009al ). At the TeV energies 
only th e upper limits are a vailable in spite of intensive moni- 
toring |Aleksic et aI.ll2010D . There are also some evidences of 
the sporadic GeV- TeV 7- ray emission from another source 
of th is type, i.e. Cyg X-l (|Albert et alj|2007l : ISabatini et all 
l20ld ). Other microquasars, such as GRS 1915+105 and 
SS433, have been also observed in the TeV energies by 
the Cherenkov telescopes but without positive result (see 
e.g. H h aroman et al.ll2005l:lHavashi et"al]|2009l :lA cero et al.l 
l2009l : ISaito et alj|2009l ; lAleksic et al.ll201ll ). The GeV- TeV 
emission has been also detected from a few other massive bi- 



nary systems, i.e. LS5039 and LSI 61 303. However it seems 
at present that they do not belong to the cl ass of micro- 
quasa rs as it was originally suspected (see e.g. iDubus et al.l 
l2006j) . 

It is expected that the 7-ray emission from microquasars 
is produced by relativistic electrons which are accelerated in 
the mildly relativistic jet. These elect rons scatter the soft 
radiation from the accretion disk (e.g. Ccrutti c t al.l 201ll), 
or th e companion star (e.g. iBednarekl |2010| : IDubus et al.l 
l20ld ). If hadrons are as well accelerated in jets of micro- 
quasars, then their 7-ray radiation should be accompanied 
by the neutrino fluxes potenti ally observable by the large 
scale neutrino telescopes (e.g. Levinson fc Waxman 200 ll : 
iDistefano et al 1 120021 : iRomero et al.ll2003l : lBednarekll2005h . 

We propose that the 7-ray emission from microquasars 
is mainly produced within the inner part of the jet in the 
anisotropic Inverse Compton e pair cascade initiated by 
relativistic electrons in the radiation field of the accretion 
disk. If the relativistic electrons are injected close to the 
disk, then the disk radiation dominates over the radiation 
from the companion star. In the example calculations we 
apply the parameters of the Cyg X-3 binary system which 
contain very hot WR type star and an object on a very com- 
pact orbit. For such luminous stars in compact binaries as 
Cyg X-3, the TeV 7-rays produced within the binary sys- 
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tern should be efficientl y absorbed in the radiation field of 
the c o mpanion star (e.g.lApparaolll984l ; IProtheroe fc Stanevl 
1 19871 : iMoskalenko et all Il993l ~ If electrons are accelerated 
close to the surface of the accretion disk, then 7-rays can 
be at first absorbed in the X-ray radiation coming from 
the disk i tself (see calculations in the case of X-ray bina- 
nes m e.g. ICarraminaiialll992l ; [Bednareklll993l; ICerutti et al 



2011 and in the case of active galaxies in Becker fc Kafatos 



1995l : ISitarek fc Bednarekl200sf ). The 7-ray spectrum emerg 



ing from the radiation of the accretion disk is modified in 
such a way that its farther absorption in the stellar radiation 
becomes negligible in most cases. Therefore, 7-ray spectra 
escaping from the compact, massive binary systems, such 
as Cyg X-3, should form in the cascade process initiated 
by relativistic electrons in the anisotropic radiation field of 
the accre tion disk as recently conside red in the case of active 
galaxies (jSitarek fc BednareHuOlOal lbT) . In this paper we cal- 
culate for the first time the angle dependent 7-ray spectra 
produced in the anisotropic IC e pair cascade initiated by 
relativistic electrons in jets of microquasars. As an example, 
we confront these spectra with the GeV 7-ray observations 
of Cyg X-3 by the Fermi-LAT telescope. 



2 A MODEL FOR GAMMA-RAY 

PRODUCTION IN IC e ± PAIR CASCADE 

We adopt the classical geometry of a microquasar in which 
a jet is launched from the inner part of the accretion disk, 
perpendicularly to its plane. The disk is formed by the mat- 
ter which is supplied by the companion star (see Fig. [T]). 
This matter is supplied to the compact object in the form 
of a fast wind. However, this wind have significant angu- 
lar momentum in the rest frame of the compact object due 
to its fast orbital motion around the massive star. There- 
fore, the formation of the relatively thin accretion disk 
in the inner part of the accretion flow is very likely. The 
standard accretion disk model is applied in which the disk 
surface temperature is described by the power law profile: 
T(r) = Ti n (r/ri I1 )~ 3 ' /4 , where Ti n is the temperature at th e 
inner radius of the disk, n n (see lShakura fc Sunvae vl ll973l ). 
For the example calculations presented below we assume 
typical values for the accretion disk, 7j n = 5 ■ 10 6 K and 
ri n = 10 7 cm. For these parameters the disk luminosity is 
equal to Ld = 4ivrf n asBT^ « 4.5 x 10 37 erg s" 1 , where <tsb 
is the Stefan-Boltzmann constant, which is typical in the 
case of X-ray binary systems. We assume that the jet has a 
conical structure with the opening semi-angle of the order 
of a = 5°. Assuming that the total flux of the magnetic field 
through the cross section of the jet is conserved, we estimate 
the magnetic field strength in the jet at the height h above 
the accretion disk on: 



B(h) = B iu rf n /(n n + hsina) 2 , 



(1) 



where _Bj n is the strength of the magnetic field at the height 
r- m above the disk. We relate the free parameter £?j n to 
the inner disk temperature by assuming some level of the 
equipartition, 77, between the energy density of the mag- 
netic field and the energy density of the isotropic black 
body radiation with the fixed temperature of Tj n : B lrl « 
4- 10 5 % /7;(ri n /10 6 K) 2 G. Note that due to the assumed tem- 
perature profile and the geometry of the radiation field close 
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Figure 1. Schematic representation of the scenario considered 
in the paper. IC pair cascades are initiated by relativistic 
electrons contained within the emission region (the blob) mov- 
ing along the jet axis. Electrons scatter in the IC process the 
anisotropic radiation coming from the accretion disk formed by 
the matter accreting from the companion star. 7-rays (the thick 
dashed arrows) are produced in the blob (the dark circle) mov- 
ing along the jet (the light shaded region) via the IC scattering 
of the thermal disk photons (the thin dashed arrow). Some of 
them escape to the observer but others either hit the accretion 
disk surface or they are absorbed in the thermal disk radiation 
and produce e+e — pairs (the thick solid line). These pairs can 
generate the second generation of 7-rays. The companion star is 
shown as the big circle on the right. The accretion disk spreads 
throughout the range of distances r- ln — r ou t and emits black body 
radiation with the temperature profile T ~ r" 3 / 4 . 



to the disk, even for r\ = 1 the magnetic field energy density 
will dominate over the energy density of the radiation field. 
For example, at the height of Inn, for r\ — 1 the magnetic 
field energy density is larger by one order of magnitude. 

In order to find out which radiation field is mostly rel- 
evant in the 7-ray production in different parts of the jet, 
we compare the energy densities of the radiation fields from 
the accretion disk and from the companion star. We assume 
typical parameters of the companion stars and of the binary 
systems as observed in two well known binaries, i.e. Cyg X-l 
and Cyg X-3 (see Fig. [2]). It is clear that the energy density 
of the accretion disk dominates up to the height comparable 
with the separation, a, of the components of the binary sys- 
tem, i.e. H ~ 10 4 - 10 5 r in ~ a. Therefore, we conclude that, 
except of the most outer parts of the jet, it is likely that 
particles accelerated in the jet will predominantly interact 
with the radiation field of the accretion disk in these two bi- 
naries. Note, moreover that the surface temperature of the 
companion star is about two orders of magnitude lower than 
the temperature at the inner disk radius. Because of that, 
we can safely neglect the effects of farther absorption and 
cascading of the 7-rays, produced in the considered model, 
in the radiation field of the companion star. 

We consider the scenario in which electrons are accel- 
erated in jets of microquasars in the classical shock acceler- 
ation mechanism. The electrons obtain the power law spec- 
trum and they are distributed isotropically in the reference 
frame of the emission region (the blob). The blob moves 
relativistically along the jet with the specific velocity Vb, 
resulting in an anisotropic distribution of the electrons in 
the disk frame of reference. The electrons interact with the 
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Figure 2. Energy density of the radiation of the accretion disk, 
with the temperature T; n = 5 ■ 10 6 K at its inner radius r; n = 
10 7 cm, as the function of the height above the disk (the solid 
curve). The energy density of the radiation from the companion 
stars: in Cyg X-l binary system (the surface temperature of the 
companion star T* = 3.2 ■ 10 4 K, the radius i?* = 1.2 ■ 10 12 cm, 
and the separation a = 2.47 -R*, the dashed curve) and Cyg X-3 
(T* = 9 • 10 4 K, R* = 2 • 10 11 cm, and a = 2.25 Ft*, the dotted 
curve) . 

radiation of the disk producing the first generation of the 
7-rays in the Inverse Compton scattering process. While the 
electrons are confined in the jet by the magnetic fields and 
move with the bulk velocity of the blob, the 7-rays produced 
by them escape from the jet into the volume above the ac- 
cretion disk. The primary 7-rays propagate in the dense ra- 
diation field of the accretion disk. Significant part of them is 
absorbed, producing the first generation of e pairs. As the 
result of the IC and 7-7 absorption processes, the IC e ± pair 
cascade devel ops in the whole volume a bove the accretion 
disk (see also lSitarek fc Bednarekll2010al lbT). In the cascade 
process considered in this paper we neglect the effects of 
bending of the trajectories of the secondary leptons by the 
local magnetic field which could be present above the disk. 
In fact, these effects can be neglected in the case of a dipole 
magnetic field above the accretion disk. Such field drops 
fast with the distance from the centre of the disk (propor- 
tionally t o oc R~ 3 ) and is important only very close to the 
disk (see lSitarek fc BednarekTl2010al ). Since the disk radia- 
tion field is highly anisotropic, the 7-ray spectra produced 
in such a cascade process strongly depend on the location of 
the observer in respect to the axis of the accretion disk. In 
the subsections below we consider the input parameters for 
such a model and show the conditions for which the 7-rays 
can both be efficiently produced and escape to the observer. 

2.1 Acceleration of electrons in the jet 

The rate of the energy gain of an electron with the energy 
E due to the acceleration process in the jet can be generally 
parametrized as: 

(dE/dt)^ = £cE/R L « 10 13 eB eV s~\ (2) 

where £ is the so-called acceleration coefficient, i?L is the 
Larmor radius of electrons, and B is the magnetic field 
strength (in Gauss) at the acceleration region. We estimate 



the acceleration coefficient on £ ~ (i'/ c ) 2 in the case of the 
acceleration process occurring in the turbulent medium of 
the jet and/or multiple shocks, where v is the velocity of 
the shock and c is the speed of light. 

The primary electrons accelerated in the jet lose their 
energy due to both the synchrotron and the IC cooling. 
Which of those dominates at a given location in the jet de- 
pends on the assumed parameters of the model, especially on 
the level of the equipartition between the magnetic field and 
the disk radiation field energy densities, rj. The synchrotron 
energy loss rate can be calculated from the formula: 

(dE/dt) Bya = -ca T U B 7e « 8.6 x KT 4 B 2 7 e 2 eV s" 1 , (3) 

where 7 C is the Lorentz factor of the electron. 

The energy loss rate due to the IC scattering is more 
difficult to compute because of the strongly anisotropic ra- 
diation field and the complicated cross section. For the con- 
sidered parameters, the electrons scatter the soft radiation 
mainly in the Klein-Nishina (KN) regime of the IC process. 
Therefore, the synchrotron energy losses can dominate even 
for relatively small values of rj. We calculate numerically the 
energy losses of electrons on these two processes as a func- 
tion of their energy and the distance from the disk surface 
(applying the full KN cross section). The energy gain and 
loss rates of the electrons are compared in Fig. [3]for different 
sets of the parameters describing the content of the jet, the 
acceleration process, and the distance from the base of the 
jet. Note that except for the case of a very weak magnetic 
field in the jet and a large distance from the base of the jet, 
the acceleration process of electrons is saturated by the syn- 
chrotron energy losses. However, just below this saturation 
energy the IC energy losses can be easily comparable or even 
much larger than the synchrotron energy losses. Therefore, 
we conclude that the significant amount of energy accumu- 
lated in primary electrons should be converted to the 7-rays 
in the IC process of the disk radiation. Closer to the base 
of the jet, where the magnetic field is stronger (see Eq. [2] 
and Eq. [3]), the energy loses grow faster than the energy 
gains, resulting in lower values of the maximum energy of 
the accelerated electrons. At H ~ r- m , electrons can reach 
maximum energies only of the order of a few to a few tens 
GeV. However, already at the distance of H ~ 10 2 r- ln , elec- 
trons can be accelerated up to ~ 10 2 — 10 3 GeV. Therefore, 
in principle, 7-rays with energies clearly above the threshold 
of even present Cherenkov telescopes should be produced in 
such a model. 

The maximum energies of the electrons, estimated 
above, are correct provided that their acceleration time is 
significantly shorter than the ballistic time of the plasma 
flow through the jet, Tb, so the electrons can be considered 
as confined at a given location in the jet. Only in such a case, 
the electrons can reach large energies close to the accretion 
disk, the source of the strong radiation. The ballistic time 
scale for the jet plasma can be estimated from, 

r b = H/rib » 3.3 x W- 4 H in //3 s, (4) 

where Hi n = H/r^ is the distance measured from the base of 
the jet in units of the inner radius of the disk, and f3 — v^/c. 
Electrons are accelerated locally within the jet (r aC c < Tb) 
for the condition, 

£acc < 3.3H in £B/p GeV. (5) 




Figure 3. Energy gains (the horizontal dashed lines), synchrotron (the dotted lines) and IC (the solid curves) losses of the accelerated 
electrons at the height H = 1 r; n (the thin lines) or H = 100 r; n (the thick lines) above the accretion disk. The equipartition parameter 
between the magnetic and the radiation fields energy density at the inner radius (see text for details) is assumed to be r\ = 1 (the left 
panels), 10 -3 (middle), and 10 — 6 (right). The calculations are shown for different acceleration efficiencies of electrons: £ = 0.1 (the 
top panels) and 0.01 (bottom). The disk temperature at the disk inner radius, r; n = 10 7 cm, is T; n = 5 ■ 10 6 K. The circles mark the 
intersections at which the energy gains balances the energy losses of electrons. 



Let us introduce some extreme parameters for this last for- 
mula. The magnetic field strength in the inner jet is of 
the order of Bi n ~ 10 4 G, even for very small rj (e.g. for 
r\ = 1CP 6 , and Hi n = 1. The typical values for the accel- 
eration parameter in the relativistic shocks are f > 0.01 
and the velocities of the microquasars jets are estimated on 
j3 = 0.3. Even for this limiting parameters, these critical en- 



ergies of accelerated electrons are of the order of E a 



10 J 



GeV. We conclude that in the model considered by us, the 
electrons can be accelerated locally in the jet even up to 
TeV energies. Note that the velocity o f the jet in Cyg X-3 
is measured in the range 0.3c — 0.8c dSpencer et al.l 19861 ; 



iMioduszewski et al.ll200ll ; iMarti et all 120011 ). While, as the 
above calculation show, the acceleration in our model can 
be considered as local, the cooling process of electrons does 
not need to occur locally in the jet. Therefore, we include in 
our numerical code the effects of the change of location of 
blob, containing the primary electrons, in the jet during its 
propagation with the jet flow. 



2.2 Absorption of gamma-rays in the disk 
radiation 

As it was discussed in Sect. [T] 7-ray photons propagating 
close to the accretion disk can be absorbed producing sec- 
ondary e ± pairs. The optical depths for 7-rays will strongly 
depend on their injection distance and the propagation an- 
gle, measured in respect to the axis of the disk (see e.g. 
ISitarek fc Bedna rck 2008; ICerutti et aljlioilh . As an exam- 
ple, we show the calculations of the optical depths as the 
function of the energy of the gamma-rays (calculated in the 
reference frame of the disk) in the case of the assumed by 



us in the Sect l2.ll accretion disk parameters, selected dis- 
tances from the base of the jet, and two injection angles of 
7-rays (see Fig. In general, the optical depths can easily 
reach values up to ~ 100. In contrary to the case of blazars, 
the typical accretion disks occurring in microquasars have 
much larger surface temperatures. This results in an effi- 
cient absorption of the 7-rays already with energies of a few 
GeV, provided that they are already produced close to the 
base of the jet, i.e. within ~ 100 n n . The optical depths also 
strongly depend on the observation angle of the accretion 
disk. These features are generally consistent with the early 
calculati ons for the accretion disks with ot her parameters 
(see e.g. ICarraminana|[l992l ; lBednare"klll99 3). As a result of 
these strong absorption effects, the 7-ray spectra escaping 
to the observer should form in the process of the IC pair 
cascade developing above the disk surface. The 7-ray spec- 
tra emerging from such a cascade should strongly depend on 
the location of the observer with respect to the axis of the 
disk, and the injection distance of primary electrons from 
the base of the jet as already consider ed in the case of the 
accret ion disks in active galactic nuclei (jSitarek fc Bednarekl 
l2010al lbl). 



3 IC e ± PAIR CASCADE SPECTRA 

In this section we investigate the IC e ± pair cascade spec- 
tra as a function of the parameters describing the model. 
We perform a calculation of the pair cascade using a M onte 
Carlo code described in (ISitarek fc Bednarekl l2010bT ). We 
expanded the original code to track also the spectra of the 
synchrotron radiation produced by electrons in the jet. The 
jet velocity is assumed to be constant in its part, where the 
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Figure 4. Optical depths of 7-rays in the radiation field of the ac- 
cretion disk with the surface temperature defined by the Shakura 
& Sunyaev profile, and the temperature T ln = 5 ■ 10 6 K at the 
inner radius of the disk rj n = 10 7 cm as a function of their ener- 
gies. The optical depths for the 7-rays propagating along the jet 
axis (9 = 0°) are marked by the thick curves and for 9 = 60° by 
the thin curves. The different curve styles correspond to differ- 
ent injection height of 7-rays above the accretion disk: H = 3ri n 
(solid), 10r; n (dashed), 30rj n (dotted), and 100 r; n (dot-dashed). 

injection of electrons occurs. The traverse size of the accel- 
eration region is assumed to be much smaller then its height 
above the disk. However we also take into account a possi- 
bility that the blob is traveling along the jet resulting in an 
elongated acceleration region (see Section \3. 2 1 . 

At first we show what is the effect of the IC e pair 
cascade process on the emerging 7-ray spectra in compari- 
son to the first generation (primary) 7-rays produced in the 
cascade and also to the case of the simple exponential ab- 
sorption law of these primary 7-rays (see Fig. [5} ■ The full 
cascade 7-ray spectra show a clear excess at energies be- 
low ~ 1 GeV (i.e. the threshold for the absorption process) 
in respect to the primary 7-ray spectra. This excess would 
have been neglected in the 7-ray spectra which are modified 
only by the simple exponential absorption law. Moreover 
the simple exponential absorption cut-off in the 7-ray spec- 
trum above ~ 1 GeV is partially filled up by the secondary 
cascade 7-rays produced by electrons with larger energies 
(compare the thick dotted and the solid curves in Fig. 5). 
The full cascade spectrum steepens sharply in comparison to 
the primary 7-ray spectrum (which corresponds to the case 
of an optically thin, complete cooling IC 7-ray production). 
The primary 7-ray spectrum has a differential spectral in- 
dex ~ 2, while the one of the full cascade spectrum above 
the break is closer to 3 — 4. The effects mentioned above 
are especially visible in the case of the large observation an- 
gles, for which the cascading effects are the strongest due 
to the favourable interaction angles between the 7-rays and 
photons coming from the accretion disk. 

3.1 Spectra for fixed distances from the disk 

We have calculated the 7-ray spectra produced in the cas- 
cade process by electrons injected in the blob moving with 
the specific velocity, Vb = 0.5c, and located at a fixed dis- 
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Figure 5. The 7-ray spectral energy distribution (SED) pro- 
duced in the IC pair cascade (the solid curves) for two ranges 
of the observation angles — 40° (the thin curves), and 60 — 75° 
(thick). The SEDs of the first generation of 7-rays produced by 
the primary electrons without any absorption taken into account 
are shown with the dashed curves. The SEDs of the 7-rays escap- 
ing from the binary system in the case of only simple exponential 
absorption law are shown with the dotted curves. It is assumed 
that the electrons are injected in a stationary blob at the height 
H = 100 n n above the disk. They have a differential power law 
spectrum with the index of —2 between 0.1 GeV and the maxi- 
mum energy obtained from the comparison of energy gains and 
losses (Sect. 2.1). The magnetic field in the jet is described by 
the equipartition coefficient r\ = 10 -6 and the decay law along 
the jet given by Eq. 1 (see text for details). Electron spectra are 
normalized to 1 erg. 

tance from the accretion disk (see Fig. [6} . In these calcula- 
tions electrons are injected into the blob at a given location 
in the jet but after the injection process the blob changes its 
distance from the disk. As the calculations of a cascade in an 
anisotropic radiation field are very time and memory con- 
suming we calculate the high energy component (from the 
IC scattering) only down to 0.1 GeV, which is the energy 
range accessible by the Fermi-LAT instrument. Presently 
there is no sensitive enough experiment to investigate the 
MeV range. Moreover in order to accurately predict the MeV 
emission (the joining point of the synchrotron and IC com- 
ponent), one would have to do additional assumption on the 
acceleration mechanism (e.g. what is the minimal energy, 
from which electrons are accelerated). 

Two scenarios for the magnetic field distribution in the 
jet are considered. In the first one, model I (defined by 
the equipartition parameter equal to 77 = 10~ 3 ), the en- 
ergy losses of electrons on the synchrotron and IC processes 
are comparable. In the second one, model II (77 = 10 -6 ), 
the energy losses on the IC process clearly dominate over 
their synchrotron losses. In the model I, the synchrotron 
spectra produced by the primary electrons in the blob al- 
most do not depend on the distance from the base of the 
disk. This feature is easy to understand keeping in mind 
that maximum energies of electrons accelerated in specific 
parts of the jet are determined also by the balance between 
the energy gains and the synchrotron energy losses. As a 
result the maximum energies of synchrotron photons do not 
change with different strength of the magnetic field. They 
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Figure 6. 7-ray SED, produced in the IC pair cascade (rightmost component, only calculated > 0.1 GeV), observed at different range 
of angles measured in respect to the jet axis: — 40° (the solid curves), 40 — 60° (dashed), 60 — 75° (dotted), and 75 — 90° (dot-dashed). 
Leftmost component of the SEDs is produced by simultaneous synchrotron radiation of the primary electrons. The primary electrons are 
injected with the differential power law spectrum with an index —2 between 0.1 GeV and the maximum energy determined by balancing 
energy gains and losses (see Sect. 2.1). They are injected from a blob moving with the velocity v\> = 0.5c at its fixed distance from the 
base of the jet: H = 1 r- ln (the left panels), 10r; n (the middle panels), 100 r; n (the right panels). The calculations have been performed 
for the acceleration coefficient equal to 5=0.01 and the equipartition parameter describing the magnetic field strength in the blob equal 
to i) = 10 — 3 (the top panels) and 10 — 6 (the bottom panels). 



depend only on the acceleration parameter. In the model II, 
the synchrotron spectra show complicated dependence on 
the distance from the disk since the maximum energies of 
accelerated electrons are determined in this case by their IC 
energy losses (see Fig. 3). The 7-ray spectra from the IC 
e pair cascade process in both models clearly extend to 
larger energies for larger injection distances from the base 
of the disk. In general, the 7-ray spectra cut-off at lower 
energies for larger observation angles due to stronger ab- 
sorption effects for 7-rays propagating at larger angles to 
the disk axis (see Fig. 4). Note the characteristic shapes of 
the cascade 7-ray spectra escaping from the disk radiation. 
They show a clear anticorrelation in energy ranges below 
and above a few GeV with the observation angle, i.e. fluxes 
above/below a few GeV are relatively large/small for small 
angles but for large angles the situation is the opposite. As 
expected, the multiwavelength spectra in such a model for 
a weakly magnetized jet (i.e. the magnetic field strength at 
the bottom of the jet B ln ~ 10 4 G) show strong dominance 
of emission at 7-ray energies. Such sources should be clearly 
detected in 7-ray energy range showing relatively weak syn- 
chrotron emission in the range from radio to X-rays. In con- 
trary, in strongly magnetized jets (Bin ~ 3 x 10 s G), the 
power emitted in the synchrotron and in the IC bumps are 
comparable. We show quantitatively that microquasars can 
show completely different spectral features depending on the 
magnetization parameter of the jet. 



3.2 Spectra from the whole jet 

We also calculate the 7-ray spectra by assuming that the 
injection of electrons occurs not only at a specific distance 
from the base of the jet but rather extends along the con- 
siderable part of the jet. This corresponds to the situation 
in which the acceleration process of the electrons operates 
through a certain period of time during which the blob 
passes specific part of the jet (see Fig. [7J. It is assumed 
that during this period the injection rate of electrons is con- 
stant but the maximum energies of freshly injected electrons 
are determined by the balance between the acceleration effi- 
ciency and the energy losses at the specific distance from the 
disk. We also follow the cooling of earlier injected electrons 
during the propagation of the blob within the jet. As an 
example, we investigate how the synchrotron and the cas- 
cade 7-ray spectra depend on the basic parameters of the 
model for two ranges of injection distances from the disk 
(H = 1 — 10 ri n and 10 — 100 n n ). As it can be expected 
from Sect. 13.11 depending on the magnetization parameter 
of the jet (77) and the acceleration efficiency of electrons (f ) 
we predict completely different shapes for the multiwave- 
length spectra at high energies. Either the synchrotron spec- 
tra dominate in the energy range from soft X-rays up to soft 
7-ray (without significant emission at ~ 100 MeV) or the 
IC 7-ray spectra dominate with relatively weak fluxes at 
the synchrotron part of the spectrum which peak usually at 
hard X-rays. Note that depending on the discussed scenario, 
the spectral index of the differential synchrotron spectrum 
changes from ~ — 2 (for the synchrotron cooling dominated 
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regime) to ~ —1.3 (for the IC cooling dominated regime) 
(see Fig. [7j). 

The increase of the acceleration parameter to the maxi- 
mum possible values expected in our model (£ ~ 0.1) allows 
the acceleration of electrons to energies with which they are 
able to produce the IC 7-ray spectra extending up to ~ 1 
TeV energies (i.e. in the range well covered by the presently 
operating Cherenkov telescopes). We suggest that some of 
the unidentified galactic TeV 7-ray sources (without obvious 
counterpart in the radio and X-rays) might be in fact low 
mass microquasars efficiently accelerating electrons. 

The IC cascade 7-ray spectra discussed above have been 
performed for the fixed velocity of the blob along the jet. In 
Fig. [8] we show how these 7-ray spectra depend on the ve- 
locity of the blob for the range of values observed in micro- 
quasars, i.e. 0.3c— 0.9c. As expected, the relativistic Doppler 
effect is important for the jets with large velocities. Then, 
the 7-ray spectra emitted at small observation angles to the 
jet axis clearly dominate over those emitted at large angles 
even in the sub-GeV energies. Therefore, in principle the 
shape of the 7-ray spectrum observed through broad energy 
range should give some independent constraints on the ve- 
locity of the 7-ray production region in microquasars. Note 
that for large velocity jets, the spectra also extend to larger 
energies, which increases the chances for a detection of such 
sources by the Cherenkov telescopes. 

3.3 Spectra for different accretion disks 

Up to now we considered only the Shakura-Sunyaev opti- 
cally thick and geometrically thin accretion disks. In fact, 
the inner part of the accretion disk can become radia- 
tive ly inefficient (the so- called advection dominated disks, 
e.g. iNaravan fc Yilll994 ). In such case, the matter in the 
inner part of the accretion disk can become optically thin 
and the temperature of the plasma in the disk can rise by 
orders of magnitude. Such disk becomes also radiatively in- 
efficient, i.e a significant part of the generated gravitational 
energy is swallowed by the black hole. Then, the rare, hot 
plasma in the inner part of the disk produces radiation field, 
which is not dense enough neither for an efficient absorption 
of 7-rays, nor for 7-ray production via the IC scattering pro- 
cess. Therefore, we neglect this inner part of the disk when 
considering the IC e ± pair cascade processes. For compari- 
son, we have calculated the optical depths for 7-ray photons, 
for different injection angles, assuming that the inner part 
of the accretion flow does not contribute to the absorption 
of the 7-rays (see Fig. Note significant differences in the 
optical depths for 7-rays both for small and for large angles 
in respect to the disk axis. The differences in the optical 
depths in respect to the previously discussed accretion disk 
are especially important at low energies due to the lack of 
the radiation from the inner, hotter part of the accretion 
disk. 

In Fig. 1101 we compare the 7-ray spectra produced 
in the IC e ± pair cascade, developing in the radiation 
field of such a trun cated disk with the case of the whole 
IShakura fc Sunvaevl disk. If the injection of the relativistic 
electrons occurs far away from the disk, the cascade 7-ray 
spectra observed for small viewing angles are very similar 
in both cases. The inner part of the accretion disk does not 
contribute significantly to the absorption process of 7-rays 



due to the small angles between the directions of the 7-rays 
and the direction to the centre of the disk. On the other 
hand, the spectra produced at a large angles are extending 
to higher energies for radiatively inefficient disk due to the 
lack of the strong absorption in the hot radiation emitted 
from the central parts of the disk. In contrary, if electrons 
are accelerated close to the surface of the disk, then the lack 
of the emission from the inner part of the disk significantly 
reduce the level of the 7-ray spectra. Also the synchrotron 
spectra of the primary electrons are modified, especially for 
the injection close to the base of the jet, due to different rel- 
ative efficiency of the two basic cooling processes considered 
in our model. Without the emission from the hot, inner part 
of the disk, the cooling of electrons dominated by the syn- 
chrotron process, is efficient also at lower energies, resulting 
in a softer synchrotron spectrum. 



4 EXAMPLE COMPARISON WITH 
OBSERVATIONS OF CYG X-3 

Cyg X-3 is the extreme case of the binary system with the 
microquasar features due to the presence of a very luminous 
WR type star (the radius i?* = 2 x 10 11 cm and the surface 
temp erature T* = 9 x 10 4 K, see ICherepashchuk fc Moffatl 
1 1994) . It is also a very com pact system (the separation of th e 
components D = 2. 257?* ICherepashchuk fc Moffatl ll994T ). 
It has been classified to belong to the microquasar class 
due to evidences of a relativistic motion. Recently, tran- 
sient 7-ray emission has b een detected from Cyg X-3 be- 
fore t he major radio flare (|Abdo et al . 2009a: Tavani et al.l 
l2009h . The emission above 100 MeV was first modeled 
with a single power law with the differentia l spectral index 
-2.7±0.05 at at±0.20 S!/s t (|Abdo et al.ll2009af ) . More recently, 
a curved powe r-law was found to describe the GeV data 
more precisely (|Abdo et al.ll201ll ). The GeV emission shows 
modulation, of the order of ~ 25%, with the maximum when 
the compact object is behind the companion star. The level 
of this modulation seems to be much smaller than those ob- 
served in the TeV 7-ray emitting bina ries such as LS 503 9 
|Abdo et al.ll2009d l and LS I +61 303 (|Abdo et al.ll2009rj ). 
Note, that the baseline emission in the Cyg X-3 light curve 
may not be well known due to the strong background from 
the Galactic disk. However, if the reported level of modula- 
tion is real, then this might suggest a different mechanism 
of the modulation of the 7-ray signal seen in Cyg X-3 from 
those ones observed in the TeV 7-ray binaries, e.g. precession 
of the disk/jet or partial contribution to the 7-ray emission 
from the outer and inner parts of the jet. The former hy- 
pothesis allows to explain the modulated and unmodulated 
7-ray signal in terms of our model alone. If the whole disk/jet 
system undergoes a precession in its movement around the 
centre of mass of the binary system, the angle at which we 
see the jet is constantly changing. As we showed above this 
will result in changing of the observed emission. Even more 
interestingly, if this is the case we should observe the vari- 
ability not only in the flux level, but also in the spectral 
shape (in particular in the cut-off energy). Note moreover 
that such an effect would be easily distinguishable from the 
variations of the spectrum due to the absorption in the ra- 
diation field of the companion star. Due to large differences 
between the temperature of the disk and of the star those 
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Figure 7. As in Fig. [6] but for the continuous injection of electrons along the jet through the range of distances from the disk: 
H = 1 — 10r; n (the top panels), 10 — 100 r; n (the bottom panels). The equipartition parameter and the acceleration efficiency are fixed 
on: rj = 1,£ = 0.01 (the far left panels) 77 = 10 _3 ,£ = 0.01 (the middle left panels). 77 = 10~ e ,£ = 0.01 (the middle right panels) and 
77 = 10 — 6 ,£ = 0.1 (the far right panels). The velocity of the blob is assumed to be 0.5c. Rightmost component of the SED is produced 
by the IC scattering in the pair cascade (only calculated > 0.1 GeV), while the leftmost component is the simultaneous synchrotron 
radiation of the primary electrons. 



spectral features would occur at completely different ener- 
gies. The latter hypothesis might be in accordance with our 
model in which acceleration of electrons occurs already close 
to the base of the jet, resulting in the 7-ray production in 
the IC pair cascade occurring in the radiation field of 
the accretion disk. In this case, the low level modulation 
of the 7-ray signal might become a consequence of the cas- 
cade developing in the radiati on field of the m assive star as 
previously considered by e.g., l|Bednarekll201(]| ). 

We apply the likely parameters of the accretion disk 
around compact object in Cyg X-3 (ri n = 2 x 10 6 cm and 
Tin = 1.5 x 10 K, consistent with the observed X-ray lu- 
minosity from the accretion disk, which is of the order of 
~ 10 38 erg s _1 ). The jet velocity i s fixed on t>b = 0.5c, within 
the range measured for Cyg X-3 (|Mioduszewski et alj|200l1 ; 
iMiller- Jones et al.|[200l ). The inclination angle of the axis of 
the binary system (and probably also of the jet) to the ob- 
server in estimated in the range i — 30° — 70° for the extreme 
param eters of the binary system (see ISzostek fc Zdziarskil 
120081 ). For these parameters we calculate the 7-ray spectra 
produced in the anisotropic IC e pair cascade as considered 
above. We compare the spectra calculated in the framework 
of our model, with the Fermi-LAT observations of Cyg X-3 
assuming that bulk of its emission comes from the inner jet 
(see Fig.fTT)). 

Investigation of the reasonable parameter range for our 
model shows that in order to reproduce the 7-ray spectrum 
observed by Fermi-LAT, the injection of electrons should oc- 
cur up to at least a few hundred ri n . As an example, we show 
the specific injection case of electrons between 1 — 300n n (see 
Fig. 11). If this condition is not fulfilled, then strong absorp- 
tion of 7-ray in the disk radiation cuts-off the GeV spectra 
at lower energies. Therefore, our mode lling lower the restric- 
tions presented bv lCerutti et al.l (|201ll ). who concluded that 
the emission region should be at least above 10 — 10 3 ri n from 
the base of the jet. The difference in our interpretation is 
due to the consideration of the full anisotropic cascade IC 



e pair cascade in our case in contrary to the simple ab- 
sorption arguments made bv lCerutti et all (|201ll ). The best 
agreement with the observational data is obtained for rather 
hard electron spectra, with the differential spectral index 
in the range 1.5 - 2, and the injection power in electrons 
equal to 2 — 8 x 10 36 erg/s. Such a value is of the order of a 
few percent of the disk luminosity. According to our model, 
the low observation angle (30 — 45°) is preferred. Interest- 
ingly, the flat part of the X-ray spectrum observed by RXTE 
above ~ 30 keV can be naturally explained in terms of our 
model as the synchrotron radiation of the primary electrons 
within the jet. Such an emission is strong enough for a rather 
large equipartition parameter (77 ~ 0.05). For low values of 
77, the hard X-ray emission requires another active mecha- 
nism which can be located e.g. in the inner disk corona (e.g. 
iHialmarsdotter et al.ll200fj| ). 



5 CONCLUSIONS 

We considered for the first time the anisotropic IC e ± pair 
cascade model for microquasars in which electrons injected 
in the inner part of the jet interact with the radiation field 
of the accretion disk. In this model, the maximum energy to 
which the electrons are accelerated in the jet is determined 
by the balance of the energy gain, and the synchrotron and 
IC energy losses, with the additional condition on the bal- 
listic escape time of the electrons being shorter then the 
acceleration time. We showed that the simplified model of a 
purely exponential absorption will result in neglecting the 7- 
ray bump at energy of ~ 1 GeV and a more steep spectrum 
in the multi-GeV energy range. As a result, the constraints 
on the location of the emission region within the jet should 
change in the full cascade scenario. We considered two possi- 
ble regimes for cooling of primary electrons, i.e. dominated 
by either synchrotron or IC process. In the framework of 
those scenarios, we investigated the expected synchrotron 
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Figure 8. As in Fig.[7]but for different velocities of the blob in the 
jet: Dj, = (emitting region is at rest, the top panels), i>(> = 0.3 c 
(the middle panels) and v b = 0.9c (the bottom panels). Injection 
occurs at the height H = 1 — 10r; n (the left panels), 10 — 100 r- ln 
(the right panels). The equipartition parameter and the accel- 
eration efficiency is fixed on: rj = 10~ 3 , £ = 0.01 respectively. 
Rightmost component of the SED is produced by the IC scatter- 
ing in the pair cascade (only calculated > 0.1 GeV), while the 
leftmost component is the simultaneous synchrotron radiation of 
the primary electrons. 



and IC spectra and their dependence on the observation an- 
gle and location of the emission region. We found out that 
for specific conditions, multi-GeV 7-ray emission can be pro- 
duced and escape absorption from the radiation of the ac- 
cretion disks in microquasars. Moreover, if the acceleration 
coefficient, £, is as large as ~ 0.1, even the TeV emission can 
be expected, providing that the jet is weakly magnetized. 
We suggest that such a scenario may explain the presence 
of some of unidentified TeV sources which do not show a 
clear radio nor X-ray counterpart. 

We also studied the effe ct of the possible tra nsition of 
the inner part of the standard lShakura fc Sunvaevl accretion 
disk into a radiatively inefficient, advection dominated disk 
on the emerging spectra. We found out that the resulting 
spectra are significantly different in two ways. If the injection 
of the relativistic electrons occurs far from the disk, the 7-ray 
spectra observed at a large inclination angle are extending 
to higher energies due to the weaker absorption. In addition, 
the 7-ray emission in the case of acceleration of the electrons 
close to the accretion disk is dimmer, due to lack of the 
strong radiation from the hot inner part of the disk. 

Finally, we also confront our model with the observa- 
tions of the Cyg X-3 microquasar performed in X-ray and 
GeV 7-ray energy ranges. We conclude that our model can 
explain the GeV 7-ray emission seen from this source, pro- 
vided that electrons are accelerated within ~ 300 ri n from 
the base of the jet with a rather hard differential spectral 



slope of —1.5 to —2. Moreover relatively small observation 
angle is required (30° — 45°). While the X-ray emission of 
Cyg X-3 is clearly dominated by the emission from the disk, 
a flat part of the SED seen by RXTE above 30 keV can be 
explained in the framework of our model as the synchrotron 
radiation of the primary electrons cooling in the jet, pro- 
vided that the equipartition parameter 77 is close to ~ 0.05. 
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Figure 9. The optical depths for 7-rays in the radiation field of the accretion disk assuming that the inner part of the disk becomes 
radiatively inefficient (the radiation from this part of the disk is neglected). The thin curves marks the results for the standard Shakura 
& Sunyaev accretion disk with temperature T; n = 5 • 10 6 K at the inner radius of the disk r; n = 10 7 cm. The thick curves are the optical 
depths for the disk in which the emission is generated only above 10r; n . Different curve styles correspond to different heights above the 
disk: H = 3r; n (solid). 10r; n (dashed), 30rj n (dotted), and 100 rj n (dot-dashed). 7-ray photons are observed at the angles: 9 = 0° (the 
left panel), 30° (middle) and 60° (right) with respect to the jet axis. 
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Figure 10. 7-ray SED produced in the IC pair cascade ob- 
served at different range of angles measured in respect to the jet 
axis: — 40° (the soli d curves), and 60 — 75 ° (dotted). The spec- 
tra from the standard IShakura fc Sunvaevl disk (as the curves in 
Fig. are shown by the thin curves, and from the truncated 
disk (without the inner part of the disk, below 10ri n ) by the 
thick curves. The primary electrons arc injected with the differ- 
ential power law spectrum with an index —2 between 0.1 GeV 
and the maximum energy determined by balancing energy gains 
and losses (see Sect. 2.1) in the range of distances from the disk: 
H = 1 - 10r in (the left panel), 10 - 100 r [n (the right panel). 
The blob is moving with the velocity v^, = 0.5c. The acceleration 
coefficient is equal to £ = 0.01 and the equipartition parameter 
is equal to rj = 10 -3 . Rightmost component of the SED is pro- 
duced by the IC scattering in the pair cascade (only calculated 
> 0.1 GeV), while the leftmost component is the simultaneous 
synchrotron radiation of the primary electrons. 
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Figure 11. SED produced in the IC pair cascade observed 
at different range of angles measured in respect to the jet axis: 
30 - 46° (solid), 46 - 59° (dashed), and 59 - 70° (dotted), for 
the disk parameters as expected from the Cyg X-3 binary sys- 
tem (r in = 2 X 10 6 cm, T in = 1.5 X 10 7 K). The disk emission 
is represented with the thin solid curve. The results of vari- 
ous observations of Cyg X-3 in the X-ray to 7-ray energies are 
sho wn with data points: B eppoSAX (the lowest energy full cir- 
cle, |^^cchi^^n[2003), INTEGRAL (the othe r two full cir- 



clesJBird et al.ll2007T ) RXTE (the empt y triangles. ISzostek et al.l 
120081) , Fermi-LAT (the empty squares. lAbdo et alj|201ll) . Elec- 
trons are injected with the differential power law spectrum with 
an index —1.75 between 0.1 GeV and the maximum energy de- 
termined by balancing energy gains and losses (see Sect. 2.1) in 
a range of distances from the disk: H = 1 — 300 r; n . The assumed 
velocity of the jet v\> = 0.5c. The acceleration coefficient £ = 0.01 
and the equipartition parameter is r\ = 0.05 (the top panel) and 
10~ li (the bottom panel). The spectra for specific ranges of the 
observation angles are normalized to the total flux seen by the 
Fermi-LAT telescope. 



